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The Casimir force between metallic plates made of realistic materials is evaluated for distances in the nanome- 
ter range. A spectrum over real frequencies is introduced and shows narrow peaks due to surface resonances 
(plasmon polaritons or phonon polaritons) that are coupled across the vacuum gap. We demonstrate that the 
Casimir force originates from the attraction (repulsion) due to the corresponding symmetric (antisymmetric) 
eigenmodes, respectively. This picture is used to derive a simple analytical estimate of the Casimir force at 
short distances. We recover the result known for Drude metals without absorption and compute the correction 
for weakly absorbing materials. 

PACS numbers: 42.50.Pq, 42.50.Lc. 73.20.Mf 



I. INTRODUCTION 

Van der Waals and Casimir forces are among the few 
macroscopic manifestations of vacuum fluctuations. Since the 
seminal paper by Casimir 1 1 ] showing the existence of an at- 
traction between two perfect conductors separated by a vac- 
uum gap, an abundant literature has been devoted to this ef- 
fect. In particular, the relevance of retardation, finite conduc- 
tivity, and finite temperature have been studied (see, e.g., 
Exhaustive lists of references can be found in several review 
papers such as i3l0|5|. 

In the last five years, the interest in Casimir forces has 
increased due to the existence of new measurements with 
improved accuracy |6, 7]. This has challenged theoreti- 
cians to quantify the corrections to the ideal case (zero tem- 
perature, perfect conductors, flat interfaces) that must be 
taken into account for an accurate comparison with experi- 
ments 18, 9, 11 [HI El El. Furthermore, the developments 
of micro-electromechanical systems (MEMS), for example, 
have shown that the Casimir effect is becoming an issue in 
nano-engineering |14, 15 1. Indeed, these short-range forces 
could seriously disturb the performances of MEMS |,16i.l . 

From a theoretical point of view, different methods exist to 
calculate Casimir forces. Casimir himself 1 1] determined the 
electromagnetic eigenfrequencies of the system and summed 
them in order to obtain the system's zero-point energy. The 
force is found by differentiation of this energy with respect 
to the geometrical distance separating the bodies I L .ItIi . In- 
genious subtraction procedures are often required to obtain 
a finite value for the Casimir energy, and realistic dispersive 
or absorbing materials can be dealt with using contour inte- 
grals over complex frequencies flsil . Another method, used 
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by Lifshitz 1 19], considers fluctuating cuiTents driven by ther- 
mal or vacuum fluctuations in the whole space. These cur- 
rents, whose spatial correlations are known through the fluc- 
tuation dissipation theorem, interact via the electromagnetic 
fields they radiate. The force is obtained by calculating the 
flux of the Maxwell stress tensor across a surface separating 
the bodies. One thus gets an integral over all possible par- 
tial wave contributions. For two parallel plates separated by 
a vacuum gap, for example, the partial waves can be labelled 
by their frequency, wave vector parallel to the interface, and 
polarization. By using clever contour deformation, Lifshitz 
greatly simplified the calculation of the Casimir force integral. 
The principal drawback of this approach is that the integrand 
can no longer be interpreted as a force spectrum. 

In this paper, we use an alternative approach and study the 
force integral over real frequencies and wave vectors. We 
show for generic materials (semiconductors and real metals) 
that in the near-field regime (separation distance small com- 
pared to the wavelengths considered), the frequency spectrum 
of the force exhibits peaks located close to surface-polariton 
frequencies. These peaks give the essential contribution to 
the Casimir force in this regime. We identify two types of 
resonant surface modes, binding and antibinding, that con- 
tribute respectively with an attractive and a repulsive term to 
the force. This substantiates early suggestions 12(1 l2lil that 
the Casimir force is due to surface modes, see also the recent 
papers by Genet et al. 1 13, 22]. 

We finally focus on materials whose dielectric function is 
modeled by a Lorentzian resonance, including a nonzero ab- 
sorption. We are able to use the qualitative suggestions men- 
tioned above and propose a quantitative estimation of the 
Casimir force in terms of coupled surface resonances. The 
dominant contribution of these resonances at nanometer dis- 
tances allows to perform exactly the integral over the mode 
frequencies, whereas the integral over the wave vector is com- 
puted to first order in the absorption. We show that the respec- 
tive contributions of binding/antibinding modes give a simple 
and accurate analytical estimate for the short-distance Casimir 
force, recovering previous results for nonabsorbing Drude ma- 
terials |Ic| l. In the corresponding Hamaker constant, we in- 
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elude eorreetions due to material losses. The accuracy of our 
results is established by comparing to numerical evaluations 
of Lifshitz theory, using tabulated data for the dielectric func- 
tions |23]. The paper concludes with a discussion of possi- 
bilities to "tune" the Casimir force that are suggested by our 
approach. 



II. SURFACE RESONANCES IN THE FREQUENCY 
SPECTRUM 

The starting point for our calculation of the Casimir force 
is Rytov's theory of fluctuating electrodynamics in absorbing 
media |24] that has been used by Lifshitz in his seminal pa- 
per 1 19|. This scheme applies to dispersive or absorbing ma- 
terials, as long as their dielectric response is linear. It has also 
been shown to provide a suitable framework for a consistent 
quantization procedure of the macroscopic Maxwell equations 
(see |25, 26] and references therein). 

In the following, we focus on the standard geometry of two 
planar half-spaces made from identical material (of local com- 
plex dielectric constant e{u!)) and separated by a vacuum gap 
of width d. In the Rytov-Lifshitz method, the Casimir force 
is computed from the expectation value of the Maxwell stress 
tensor at an arbitary position in the gap. At zero temperature 
and after subtraction of divergent contributions, Lifshitz gets 
a force per unit area given by 1.1 9il 



F = 



/ (U.LO) 



2tt 



F{u,uj) = _^Iml f 2^ - 

\ A' = s, p 



r2(u,w)e-2'^''^/'= 




where v = {u^ — 1)^/-^ (Imt; < 0), and is the Fresnel re- 
flection coefficient for a plane wave with polarization ^ and 
wavevector K = oju/c parallel to the vacuum-medium in- 
terface. We use the convention that an attractive force cor- 
responds to F > 0. We note that Rytov's approach allows 
for an easy generalization to different media held at different 
nonzero temperatures. The radiation force on a small polariz- 
able sphere above a heated surface has been discussed previ- 
ously in iEtIi . Results for the non-equilibrium Casimir force 
will be reported elsewhere. 

Lifshitz evaluated the integrals Q by deforming integra- 
tion contour in the complex plane to arrive at an integral 
over imaginary frequencies u; = i^. The integration then re- 
quires the continuation of the dielectric function from real- 
frequency data to e(i^), using analyticity properties as dis- 
cussed in |9, 10]. We follow here a different route and con- 
tinue to work with real ui and u, taking advantage of the fact 
that Lifshitz' results provides us with an expression for the 
frequency spectrum F{uj) — J F{u, uj)du/2TT of the Casimir 
force. Note that the force spectrum is more difficult to define 
in a calculation based on mode summation, see, e.g., ]28, 29]. 

For a polar material like SiC, the spectrum of the force is 
dominated by narrow peaks in the UV and in the IR (Fig.Q 
when the distance d is reduced to the nanometer range. These 
peaks can be ascribed to the surface phonon polaritons in the 
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FIG. 1 : Contributions of s and p polarized, propagating and evanes- 
cent modes to the force spectrum (Eq. J2j integrated over the 
wavevector it). Distance d = 10 nm. Material: SiC, dielectric func- 
tion taken from tabulated data 1 23]. The corresponding surface reso- 
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FIG. 2: Contributions of s and p polarized, propagating and evanes- 
cent modes to the force spectrum (Eq. J2j integrated over the 
wavevector u). Distance d = 10 nm. Material: aluminum, described 
by tabulated optical data 123.1 . 



IR and to surface plasmon polaritons (SPP) in the UV. The 
largest contribution comes from the UV surface plasmon po- 
lariton even though larger losses make it broader The large 
difference between the UV and the IR contributions in Fig.^ 
is due to the factor u)'^ in Eq. (|2jl. In Fig.|2j we plot the spec- 
trum of the force between two aluminum half-spaces, using 
tabulated data for the dielectric function ]23]. The dominant 
contribution to the force is clearly due to the surface plasmon 
polaritons. Indeed, the frequency of the peaks corresponds to 
the frequency Vl of the asymptote of the SPP dispersion rela- 
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FIG. 3: Dispersion relation [Eq. J3}] of the surface plasmon polariton 
on a flat interface vacuum/aluminum. The dielectric function is taken 
from the data tabulated in Ii23i1 . We plot the real part of uj versus the 
real part of the parallel wavevector K — uu> / c. 
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where the sign of the square root is chosen such that 
Rewspp > 1. It is seen in Eq. (|3} that the frequency fi is 
given by the condition Ree(ri) = — 1. This corresponds to a 
large increase of the density of states and therefore to a peak 
in the energy density l3ll l32|. The polarization dependence 
of the force spectrum provides a second argument in favor of 
a surface plasmon polariton. In Figs. [2 121 we have separated 
the contributions to the spectrum according to the mode polar- 
ization (s or p). The modes in the cavity can further be classi- 
fied into evanescent (surface) modes (u > 1) and propagating 
(guided) modes (0 < w < 1). Among the four contributions 
it is seen that the leading one comes from the p-polarized sur- 
face modes, of which the SPP is a special case. 

It is worthwhile pointing out that for a perfectly conducting 
metal, the spectrum of the force would be completely different 
because of the lack of SPP. The usual picture of the Casimir 
effect in that case is based on the modified density of states 
for propagating waves between the two plates. This picture 
includes only what we have called guided modes and ignores 
surface (or evanescent) modes. 

We observe from Figs lll2l that the contribution of the force 
is either positive or negative depending on the frequency. We 
analyze this behaviour in the next section. 



III. BINDING AND ANTIBINDING RESONANCES 

In order to further analyze the role of SPPs for the Casimir 
force, we plot in Fig. |^ the integrand F{u, uj) as given by 
Eq. (13 for two aluminum half-spaces separated by a distance 
of d = 10 nm. Two branches emerge with dominant contribu- 
tions, the higher-frequency branch yielding a negative contri- 
bution whereas the lower branch gives a positive (attractive) 
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FIG. 4: (a) Wavevector resolved spectrum of the Casimir force 
[Eq. ^2)] in the (u, uj) plane between two aluminum half spaces sep- 
arated by a distance of 10 nm. The frequency of the flat asymp- 
tote corresponds to the peaks of the force spectrum Fig. |2| Light 
(dark) areas: attractive (repulsive) force, (b) Resonant denominator 
1/|1 — Tp e"'^""'*^'^!^ in the (u,lu') plane, the grayscale giving the 



logarithm to base 10. The dispersion relation of the coupled surface 
resonance corresponds to the light areas; dark area: dispersion re- 
lation for a single interface [Eq. J3j]. Dielectric function extracted 
from tabulated data 1 23]. The inset sketches the magnetic field of the 
coupled surface resonances (antisymmetric and symmetric combina- 
tions). 



contribution. These two branches are reminiscent of the dis- 
persion relation of a SPP on a two interfaces system. It is 
given by the complex poles of the reflection factor of the two 
interfaces system in the (u, uj) plane: 



1 



^2 ^-2ujvd/c 







(4) 



In order to illustrate the influence of the SPP dispersion re- 
lation on the force, we plot in Fig. EJ) the quantity 1/|1 — 
j-^Q-'^^'''d/c^2 jjjg j.g^j (y^^) plane. Comparing Figs.|4j3 

and|^, it is clearly seen that the main contribution to the force 
is due to the SPP. In addition, we observe on Fig.|4j) a dark line 
which corresponds to minima of 1/|1 — Tp e~2c^ud/c|2 
can be attributed to very large values of the reflection factor 
Tp. Thus, the dark line is the dispersion relation of the SPP on 
a single flat interface. Note that the Casimir force shows no 
prominent feature in this region. 

In Fig.|5j we plot the force for a spacing d — 100 nm: the 
two branches tend to merge with the flat interface dispersion 
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FIG. 5: Same as FigO but for a separation d — 100 nm. 



relation. The following interpretation thus emerges: when the 
surfaces approach each other, the overlapping of the two SPP 
leads to a splitting of the polariton frequencies |33, 34]. The 
frequency splitting can be found from the solutions of Eq. (0} 
which are implicitly defined by (see also L22.I ) 



rp{u,uj) = ±e' 



jvd/c 



(5) 



The signs correspond to either symmetric or antisymmetric 
mode functions (for the magnetic field), as shown in Ap- 
pendix A and sketched in Fig. |4j). The symmetric (antisym- 
metric) branch corresponds to a lower (higher) resonance fre- 
quency, respectively, similar to molecular orbitals and tunnel- 
ing doublets |35|. These branches contribute with opposite 
signs to the Casimir force, due to the identity 



2rp(tj, u) e 



-2uvd 



1 — rp(w, u) e~ 



jvd 



1 + rp(ijj, u) e 



—uvd ' 



(6) 



where the first (second) term is peaked at the symmetric (an- 
tisymmetric) cavity mode. The symmetry of the resonance 
mode function hence determines the attractive or repulsive 
character of its contribution to the Casimir force. We show in 
Appendix A by evaluating explicitly the Maxwell stress ten- 
sor, that symmetric modes are binding as in molecular physics. 

We note that the splitting in Eq. (|6} of the force spectrum 
gives meaningful results also after integration because for 



evanescent waves, both terms converge separately. We also 
point out that for a complex permittivity e(cj) (as required by 
the Kramers-Kronig relations for a dispersive material), the 
SPP dispersion relation necessarily moves into the complex 
plane and is never satisfied in the real (m, ijj)-plane, thus ex- 
cluding any singularities of the integral ([T). 



IV. SHORT-DISTANCE LIMIT 

The short-distance behaviour of the Casimir force between 
non-perfect metals has been computed in ^ fioll using tab- 
ulated data for the dielectric function and integrating over 
imaginary frequencies. We show here that these results can 
also be recovered with a real frequency calculation. In par- 
ticular, we prove that the interaction between SPPs across 
the vacuum gap quantitatively accounts for the short-distance 
Casimir force derived in thus completing qualitative dis- 
cussions put forward by Gerlach [21] and Genet, Lambrecht, 
and Reynaud Ell . 

For definiteness, let us adopt a Lorentz-Drude model for the 
dielectric function 



1 



(7) 



with resonance frequency cjq and damping coefficient 7. The 
corresponding plasma frequency is [2(il^ — o^q]^/^. With this 
convention, the large u asymptote of the SPP dispersion Q 
occurs at w « f2. This model can be used to describe either 
dielectrics or metals when loq — 0. In the region of large 
wavevectors, the p-polarized reflection coefficient has a pole 
at ri: 



u > 1 : rp(w, u) 



eiu) - 1 



02 



■^0 



(8) 



e(cj) + 1 - 170; 

From FigsQIJ] we know that the force is significant only in a 
range around the SPP resonance. It follows that the model for 
e{uj) is needed only in this limited range. We have checked 
that Eq.|8]with o^o = is well suited to describe the reflection 
data computed from tabulated data for aluminum. Note that 
the results of the fitted parameters {fl and 7 are indicated in the 
caption of fig|3 differ from the usual bulk plasma frequency 
and damping rates that we would get from a fit over the entire 
spectrum. 

We have checked that this formula is well suited to describe 
the reflection coefficient computed from tabulated optical data 
in the frequency region around the SPP resonance. For alu- 
minum, we get a good agreement with the values given in the 
caption of Fig|6l These values do not correspond, of course, to 
the usual bulk plasma frequency and damping rates that enter 
in the Drude model of the dielectric function at low frequen- 
cies. 

With this form of the reflection coefficient, Eq. (|5jl yields 
the following dispersion relation for the (anti)symmetric SPP 
resonances, neglecting for the moment the damping coeffi- 
cient 7: 



(9) 
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We have used v ^ u for 1. For large u, we solve by itera- 
tion and find that uj± ^ fl. As announced above, the symmet- 
ric mode (upper sign) occurs at a lower resonance frequency. 

To derive an analytical estimate for the Casimir force, we 
retain in Eq. (|2ji only the contribution of p-polarized, evanes- 
cent waves, containing the SPP resonance. Introducing the 
new variable x = uovd/c, we get using the identity (|6|l 



F 



Im / duj I x^dx X 



Jo 



V 



(10) 



where A = ±1 corresponds to symmetric (antisymmetric) 
modes, respectively. The integral is dominated by the range 
a; ~ 1 and ut ^ il. To leading order in ild/c — > 0, we can 
thus use the asymptotic form of Vp valid for large u given by 
Eq. (|8j. Performing the integral over uj analytically and in- 
cluding damping to first order in 7/^7 yields 



F = 



E 

A = ±l 



dxx X 



Xze 



7AZ e ^ 



2^1 - Aze- 



27rrj(l - Azc-^) 



(11) 
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FIG. 6: Comparison of different expressions for the Casimir force 
between aluminum surfaces. We plot the ratio F{d) / Fcas{d) where 
-Fcas(rf) = ?!.c7r^/(240d^) is the Casimir force for perfect mirrors. 
Solid line: numerical integration of Eq. Q, using tabulated optical 
data Short-dashed line with circles: same with a model 

dielectric function of Drude form [Eq. 0] with luq — 0, — 
1.66 X 10^^ s~^, and 7/0 — 0.036. These parameters have been ob- 
tained from a plot of the reflection coefficient (e{uj) — 1)/ {s{uj) + 1) 
based on the tabulated data that has been fitted to the form given in 
Eq. Jsj. Long-dashed line: short-distance asymptotics il2\ with the 
same values for uo, ^l, 7. 



where z — 1~ Wq/JI'^. This result shows clearly that symmet- 
ric and antisymmetric modes give Casimir forces of opposite 
sign. The first term in the parenthesis can be computed by ex- 
panding the square root in a pow er series in Xz e^^, leading 
to an infinite series given in llOl l22ll . The second term, the 
correction due to damping, can be integrated in terms of the 
polylogarithmic function, so that we finally have 



F 



a{z) 



where 



n=l 



lU,{z^) 



. (4n-3)!! 
n3(4n- 2)!! 



and 



(12) 



(13) 



(14) 



For completeness, we give the asymptotic series for small 

LDq/^I (z 1) 

a{z) w 0.1388 - 0.32 + 0.4 (15) 



03(2^) « C(3)--(l-z) 



(16) 



3- — -21og(2(l-z)) 



{l~zf 



with C(3) ~ 1.202. (The coefficient of the second order term 
in Eq. (I15> is only accurate up to a logarithmic correction.) 



Our result Eq. (I12> for the short-distance Casimir force 
agrees with the formula given in fioi l22il in the special case 
7 = 0, Wo = (lossless Drude model). A very similar expres- 
sion has been found in |^^. We compare Eq. \\2\ in Fig. |6] 
to the full integral Eq. (|3 for the case of aluminum: it turns 
out to be quite accurate for distances d < 0.1 Aspp where 
Aspp = 115 nm is the wavelength of the SPP with the largest 
frequency |36]. In the case of aluminum, the first order cor- 
rection in 7/r2 is 2.5% of the zeroth order value of the force. 
The plot also shows that for the numerical integration, the tab- 
ulated data and the Lorentz-Drude model with parameters 
fitted around the surface resonance give very close results over 
a large range of distances. This is another indication that the 
short-range Casimir force between real metals is dominated 
by a narrow frequency range. Differences of the order of a 
few percent appear at large distances where the Casimir force 
is dominated by the low-frequency behaviour of the reflection 
coefficient that is not accurately modelled with the fitted pa- 
rameters. 

We finally note that the correction of order ^ /Vl derived 
here introduces the effects of losses and must not be con- 
fused with the correction due to a finite real permittivity 
This is already taken into account by the finite value of the 
plasma frequency Q, and is responsible for the emergence of 
the short-distance regime d ^ Aspp where the Casimir force 
~ 1/d^ \ \9\. At large distances, a finite f2 leads to a small 
correction to the well-known Casimir force ^ l/d^ between 
perfect conductors L2..9. .10.1 . 
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V. CONCLUSION 

We have pointed out that the Casimir attraction between 
reahstic materials can be quantitatively understood, at short 
distances, in terms of the interaction between electromagnetic 
surface plasmon (or phonon) polaritons. The modes overlap 
across the vacuum gap and split into symmetric and antisym- 
metric combinations which contribute with different signs to 
the Maxwell stress tensor and hence to the Casimir force. 
We discussed in particular the short-distance regime of the 
Casimir force where F — H/ <fi and have given an analytical 
formula for the Hamaker constant H. We recover previous re- 
sults for nonabsorbing materials and evaluate a correction due 
to absorption. Our results have been validated by comparing 
to a numerical calculation based on Lifshitz theory. 

The approach presented here has the advantage of linking 
in a transparent way the Casimir force to the actual physi- 
cal properties of the material surface. This suggests the pos- 
sibility of engineering the surface plasmon polariton disper- 
sion relation to modify the Casimir force. Indeed, as it has 
been shown, the Casimir force at short distances is entirely 
due to the interaction between surface polaritons. Magnetic 
materials which exhibit Casimir repulsion ll37i| an d support 
s-polarized surface waves when Re /i < —1 |38l| are good 
candidates. The folding of the dispersion relation in recip- 
rocal space by a grating, known to change the surface wave 
behaviour 1 39] could also lead to a substantial modification of 
the Casimir force. 

Acknowledgments. — This work has been supported by the 
bilateral French-German programme "Procope" under project 
numbers 03199RH and D/0031079. 



APPENDIX A: ANGULAR SPECTRUM ANALYSIS 

In this appendix, we compute the Casimir force in terms 
of an angular spectrum representation of the electromagnetic 
fields that is adapted to the planar geometry at hand. 

Letting the vacuum gap occupy the region —d < z < 0, we 



can expand the electric field in the gap as 

fj. = s,p 

i;^(K)e''=^(^+'')e+) e'^'^ (Al) 

where K = {k^, ky) is the component of the wavevector par- 
allel to the interfaces and k^ = \/ {uj/cY — K"^ its perpendic- 
ular component. The (p = s, p) are unit polarization vec- 
tors, and Et^CK) are the amplitudes of up- and downwards 
propagating plane waves. A similar expansion holds for the 
magnetic field H(x, us) with amplitudes iJj!(K). We get the 
averaged Maxwell stress tensor by integrating incoherently 
over the contributions r^^(K) of individual modes. For the 
particular case of a p-polarized evanescent mode (K > u), 
we get by straightforward algebra 

rP,(K) = 2^iov' Re [i/P*(K)ffP(K)] . (A2) 
The up- and downward propagating amplitudes are of course 
related via the reflection coefficient from the upper interface. 
Taking the phase references in Eq. iAl\ into account, we have 

ijP = rp e^'^^'^Hl = rp e~""''/"iJP w ±iJP , (A3) 

where the last equality applies in the vicinity of the cou- 
pled surface resonances defined by Eq. (|5jl. The condition 
Tp e^"'"'/'^ = +1 thus corresponds to a symmetric magnetic 
field distribution on both interfaces, because = . In 
addition, with our sign convention, this mode gives an at- 
tractive contribution proportional to +|iJP p to the stress ten- 
sor (IA2> . The opposite is true for antisymmetric modes. 

The sign of the Casimir force due to the coupled polariton 
modes can also be understood in terms of the charge densities 
excited on the surfaces, as pointed out by Gerlach [21]. These 
can be found from the normal component of the electric field. 
For a symmetric mode, we get surface charges with opposite 
sign, hence an attractive force, while an antisymmetric mode 
corresponds to equal surface charges. 
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